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ABSTRACT 


For  heavy  nuclei,  very  few  measurements  of  the  nuclear  giant  multipole 
resonances  have  been  conducted  using  inelastic  electron  scattering.   In 

this  work,  measurements  were  made  to  identify  the  giant  multipole  resonances 

197 
of    Au.   Five  inelastic  scattering  spectra  were  obtained  and  analyzed  for 

beam  energies  and  angles  of  75  MeV  /  75  ,  75  MeV  /  90  ,  75  MeV  /  120  , 

90  MeV  /  75  ,  and  90  MeV  /  90  .   In  each  spectrum,  giant  multipole  resonances 

were  observed  at  excitation  energies  of  7.1,  9.1,  10.9,  and  14  MeV.   The 

strengths  and  multipolarities  were  found  to  be  B(M2)  =  15,000e  fm  , 

B(E2)  =  8,000e2fm4,  B(E2)  =  9,000e2fm4,  and  B(E1)  =  90e2fm2,  respectively. 

The  full  widths  at  half  maximum  for  these  resonances  were  found  to  be, 

respectively,  3.1,  2.9,  3.1,  and  4.9  MeV. 

The  resonances  found  at  10.9  and  14  MeV  confirm  the  only  other  measure- 
ments made  on  gold  |  4 1 .   The  resonance  found  at  7.1  MeV  agrees  with  those 
findings  concerning  the  magnetic  character  of  the  transition  but  differs  in 
the  multipolarity  assignment;  however,   our  assignment  may  also  be  consid- 
ered tentative  since  insufficient  data  were  obtained  at  the  backward  scat- 
tering angles.   The  resonance  found  at  9.1  MeV  has  not  been  previously 
reported. 
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I.   INTRODUCTION 

The  giant dipole  (El)  resonance  has  been  well  known  for  some  time. 
Recent  experiments  using  inelastic  electron  scattering  have  revealed  the 
existence  of  other  giant  resonances  in  a  number  of  nuclei  and  in  this  paper 
we  report  our  investigations  about  these  resonances. 

Many  inelastic  scattering  measurements  have  been  made  in  the  past  15 
years;  however,  the,  giant  resonances  considered  in  this  paper  appear  as  small 
peaks  superimposed  on  a  large  continuous  spectrum  of  electrons  called  the 
elastic  radiation  tail.   Low  background  and  high  accelerator  beam  stabil- 
ity are  required  to  detect  these  resonances  and  it  is  only  recently  that 
they  have  been  investigated.   Studies  using  Ce,  La,  Pr  [  Refs.  1,2]  and 
Pb   Ref.  3   indicate  that  heavy  nuclei  tend  to  have  an  El  resonance  at 
about  13  to  15  MeV  and  an  E2  resonance  at  about  10  to  12  MeV.   Other  reso- 
nances, some  of  which  were  magnetic  in  character,  were  observed  at  excita- 
tion energies  ranging  from  7  to  9  MeV.   Unpublished  work  at  Darmstadt 
[Ref.  4 J  on    Au  has  also  confirmed  this  trend:   an  El  resonance  was  found 
at  10.8  MeV;  and  one  other  resonance  was  located  at  7.1  MeV  and  tentatively 
identified  as  Ml.   Another  recent  experiment  on   Zr  |Ref.  5  ,  found  a  giant 
resonance  (E0  or  E2)  at  28  MeV  excitation  in  addition  to  the  usual  E2  and 
El  transitions  near  15  MeV. 

The  investigations  reported  here  were  the  first  attempt  to  measure 
giant  resonances  at  the  NPS  LINAC.   The  procedures  and  data  handling  tech- 
niques necessary  for  the  study  of  giant  multipole  resonances  by  inelastic 
electron  scatterings  have  been  developed  and  refined. 

Experimental  values  for  the  inelastic  scattering  form  factors  were 
obtained  using  incident  energies  of  75  and  90  MeV  and  scattering  angles  of 
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75,  90,  120,  and  150  degrees.   The  results  of  this  work  agree  with  other 
reported  works  |Ref.  4   in  that  the  previously  reported  resonances  were 
observed.   A  previously  unreported  resonance  was  found  at  9.1  MeV  with 
B(E2)  =  8,000e  fm  .   Additionally,  preliminary  indications  for  the  exis- 
tence  of  an  E2  resonance  in  the  vicinity  of  23  MeV  were  obtained. 

This  thesis  is  organized  into  five  sections  and  one  appendix.   In 
Section  II  a  brief  review  of  inelastic  scattering  and  giant  resonance 
theory  is  presented.   Section  III  briefly  describes  the  NPS  LINAC  and  the 
experimental  procedure  used  by  the  authors.   Section  IV  discusses  the  meth- 
ods employed  for  data  reduction.   Topics  include  radiation  corrections, 
calculations  of  the  elastic  radiation  tails,  procedures  used  in  the  loca- 
tion and  characterization  of  the  giant  resonances,  and  error  analysis. 
Experimental  results  and  conclusions  are  presented  in  Section  V.   The  cal- 
culations involved  in  this  work  are  quite  complex  and  cannot  reasonably  be 
performed  without  the  use  of  various  sohpisticated  computer  programs.   The 
appendix  contains  the  program  listings  and  detailed  explanations  concerning 
their  use. 


II.   THEORY 

A.   INELASTIC  ELECTRON  SCATTERING 

When  a  relativistic,  mono-energetic  electron  beam  is  incident  on  a 
thin  target  most  electrons  pass  through,  but  some  collide  elastically  with 
target  nuclei  and  a  few  collide  inelastically,  thereby  exciting  nuclei  to 
higher  allowed  energy  states. 

The  energy  distribution  of  the  scattered  electrons  is  determined  ex- 
perimentally at  various  scattering  angles.   Most  of  the  electrons  observed 
at  a  given  energy,  E  ,  less  than  the  incident  energy,  E. ,  are  due  to  elas- 
tic collisions  with  nuclei  accompanied  by  energy  loss  in  passing  through 
the  target.   Energy  loss  mechanisms  include:   (1)  ionization  losses, 
(2)  emission  of  a  photon  during  the  scattering  event,  and  (3)  bremsstrah- 
lung  while  passing  through  the  target  both  before  and  after  the  elastic 
scattering  event.   This  "elastic  tail,"  and  any  other  "background"  affect- 
ing the  observed  spectrum,  must  be  calculated  and  subtracted  from  the  ex- 
perimental energy  spectrum  to  obtain  the  inelastic  spectrum.   The  resulting 
unfolded  inelastic  spectrum  consists  of  the  inelastic  peaks  and  their 
associated  tails.   Each  inelastic  peak  corresponds  to  an  excited  state  of 
the  nucleus. 

Inelastic  collisions  are  classified  into  three  types,  according  to  the 
mechanism  of  interaction  with  the  nucleus: 

1 .   Surface  Nucleon  Excitation 

One  or  more  surface  nucleons  are  excited  with  excitation  energies 
ranging  from  several  KeV  to  a  few  MeV.   The  higher  energies  tend  to  be 
associated  with  collective  excitations  and  have  larger  characteristic 
resonance  widths. 
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2.  Giant  Resonances 

Many  nucleons  participate  in  oscillations  causing  the  giant  reso- 
nances.  These  resonances  have  very  large  natural  widths  and  are  found 
typically  in  the  region  of  ten  to  20  MeV  excitation  for  heavy  nuclei.   This 
large  resonance  width  occurs  because  the  excitation  energy  is  large  enough 
that  the  excited  nucleus  can  decay  by  particle  emission.   The  short  life- 
times associated  with  particle  emission  result  in  the  large  widths  because 
of  the  Heisenberg  uncertainty  relationship  between  energy  and  time. 

3.  Quasi-Elastic  Resonances 

In  this  process  a  nucleon  is  ejected  from  the  nucleus  during  the 
collision.   The  incident  electron  interacts  with  essentially  only  one 
nucleon.   These  "quasi-elastic"  collisions  produce  a  very  wide  peak  in  the 
spectrum  at  an  excitation  energy  of  a  few  hundred  MeV  for  heavy  nuclei. 

B.   GIANT  RESONANCES 

According  to  the  Generalized  Goldhaber-Teller  Model  I  Ref .  6l ,  nuclei 
are  capable  of  "collective  oscillations"  in  which  all  the  nucleons  oscil- 
late inside  the  nucleus.   There  are  four  separate  fluids  in  motion:   neutrons 
with  spin  up,  protons  with  spin  up,  neutrons  with  spin  down,  and  protons 
with  spin  down.   Each  fluid  is  characterized  by  the  spin  and  the  isospin 
of  its  particles. 

The  motion  of  the  four  fluids,  two  of  which  are  180   out  of  phase  with 
the  other  two,  can  be  characterized  by  two  numbers,  S  and  T.   If  protons 
and  neutrons  move  in  phase,  T=0;  if  they  move  180   out  of  phase,  T=l.   If 
particles  with  spin  up  are  in  phase  with  those  with  spin  down,  S=0;  if  they 
are  180  out  of  phase,  S=l.   A  transition  to  an  S=l  state  changes  the  mag- 
netization and  is  therefore  called  a  magnetic  transition.   A  transition  to 
an  S=0  state  is  called  an  electric  transition. 
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The  electric  and  magnetic  fields  of  nuclear  excited  states  are  described 
in  terms  of  multipole  expansions.   A  transition  is  labeled  according  to  the 
term  in  the  multipole  expansion  which  is  predominantly  affected.   For  ex- 
ample, a  transition  from  the  ground  state  which  increases  the  magnetic  (or 
electric)  dipole  term  in  the  electromagnetic  multipole  expansion  is  labeled 
Ml  (or  El).   If  the  electric  (or  magnetic)  quadrupole  term  is  involved,  the 
transition  is  labeled  E2  (or  M2).   An  electric  monopole  transition  is  labeled 
EO.   There  are  no  magnetic  monopole  transitions.   In  general,  L  units  of 
angular  momentum  are  required  for  an  EL  or  ML  transition  from  the  ground 
state. 

C.   NUCLEAR  FORM  FACTORS 

1.   Development  of  the  Form  Factor 

In  the  plane  wave  solution  to  the  non-relativistic  wave  equation 
we  find  that  the  differential  elastic  scattering  cross  section  can  be  ex- 
pressed as  the  magnitude  of  the  square  of  some  function  of  scattering  angle, 
f  (.6)  ,  where 

fW)  =  -4V/U(r)eU,?d3r, 

-v    2m   __  _ 

and  U(r)  =  ^2  V(r),  where  V(r)  is  the  potential  of  the  electron  in  the 

field  of  the  nucleus,  and  k  is  the  wave  vector  of  the  incident  wave. 

Replacing  U(r)  by  the  related  electrostatic  potential,  <p(r) ,    through  the 

relation  V(r)  =  q  0(r) ,  we  have 


fM=+flf/ei*'r,?H?)dV- 


We  can  relate  0(r)  to  the  charge  distribution  p(r)  through  Poisson's  equa- 
tion; taking  the  fourier  transforms  and  equating  integrands  leads  to  the 
result 
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f(0)=  ■W(«D,-£L'<*>. 


with 


i>(-)=  (iV)S>(f)e       dV  =  UV)'ZeF(s). 


Hence , 


We  can  now  solve  for  the  differential  cross  section  to  obtain 

For  a  point  charge  it  can  be  shown  that  -gj-r  ^7^y~  »  so  we  have 

do-ld£.)    |F(k)|2 

dn    ida/pt1      '  • 

We  call  F(k)  the  form  factor  for  the  charge  distribution  with 

F(*)=  i/pfie^'tff.  (ii-D 


For  inelastic  scattering,  the  charge  distribution  changes  during  the 
transition.   In  this  case,   p(r)  is  replaced  by  the  four-vector  Dirac  cur- 
rent operator,  J„ ,  associated  with  the  transition  involved. 
2.   Mott  Cross  Section 

The  differential  cross  section  per  unit  solid  angle  for  elastic 
scattering  of  a  relativistic  electron  by  a  point  nucleus  of  charge  Z  was 
shown  by  Mott  [Ref.  7]  to  be: 

Ida] -    Z2e4COS2fl/2   p  =  Ida]  n 

[dnj^    4E?sin40/2     K     IdnJmott'      ' 

where  R  is  the  recoil  factor  defined  by  the  relation 
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RS     1+__ 2Ei 

1     Mc2sin26?/2 


Using  the  concept  of  the  form  factor  we  are  able  to  generalize  this  result 
to  finite  charge  distributions  to  obtain 


Id  a\         Id  a 


|k*)|-r. 


Ida/      \dn/mott 

For  heavy  nuclei  where  E.  <  100  MeV,  R  =  l;  therefore 


lda\        lda\  |F(^|  2) 


Ida/     |da/  mott 

If  this  expression  is  used  with  the  form  factor  given  by  equation  II-l, 
we  have  the  Born  approximation  for  elastic  scattering.   For  large  Z  and  low 
energy  this  approximation  is  inaccurate  and  the  cross  section  must  be  cal- 
culated by  partial  wave  analysis  [Ref.  8J  . 

3.   Electron  Scattering  from  a  Proton 

The  electromagnetic  interaction  Hamiltonian  between  the  electron 
and  a  point  proton  is  expressed  in  four-vector  notation  as  H^i=-J^[A  . 
Calculating  JL  in  terms  of  the  Dirac  theory  operators   [Ref.  9]  gives 

H;,=  IHe)«>i)X)(-ieS)«OT,). 

Solving  the  Dirac  wave  equation  yields  the  point  Dirac  proton  solution: 

jck=  Ids.]  .R.(1  +  Jl^an'0/2)  =  Ids]  -R-  F,(q2)f 

da      idajmott  p  Ida/ mott 

The  first  term  in  the  parentheses  corresponds  to  the  electric  interaction 
and  the  remaining  term  corresponds  to  the  magnetic  interaction.   General- 
izing  to  finite  charge  distributions  and  adding  another  term,  F2 (q  ),  to 
account  for  the  anomalous  magnetic  moment  distribution  we  obtain 


14 


<kJ<k\        .r.HfI'+-9 


q* 


3Hi^)mott:R-^l  +^*>  +^tan^|F1+„F,|  ). 

The  first  two  terms  account  for  the  electric  interactions  and  the  last  term 
accounts  for  the  magnetic  scattering;  q2  is  the  square  of  the  four-momentum 
transfer  and /z  is  the  proton  magnetic  moment  in  nuclear  magnetons. 

This  result,  although  it  does  not  apply  directly  to  these  investigations, 
does  show  the  presence  of  two  form  factors,  due  to  scattering  from  the  charge 
and  magnetic  moment  of  the  proton. 

4.   Plane  Wave  Born  Approximation  (PWBA) 

The  Born  approximation  is  used  to  calculate  inelastic  scattering 
cross  sections  based  on  an  assumed  nuclear  transition  charge  distribution. 
The  PWBA  uses  first  order  perturbation  theory  to  calculate  cross  sections 
for  interactions  between  an  electron  and  the  electromagnetic  field  of  the 
nucleus.   A  multipole  expansion  is  used  to  describe  the  fields.   The  tar- 
get nuclei  are  assumed  to  recoil  non-relativistically.   The  PWBA  assumes 
plane  wave  solutions  for  the  incident  and  scattered  electron. 

For  Ei5  Ef»Ex,  Isabelle  [Ref.  id]  has  shown  that  the  differential 
inelastic  cross  section  per  unit  solid  angle  is  given  by 


idft/£L    \dnj  moft 


and 


where  FT  and  F,  are  the  form  factors  for  longitudinal  and  transverse  photon 
absorption.   For  elastic  scattering,  the  PWBA  results  are  equivalent  to 
equations  11-1  and  II-2. 
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5.   Distorted  Wave  Born  Approximation  (DWBA) 

For  heavy  nuclei,  the  coulomb  field  surrounding  the  nucleus  dis- 
torts the  electron  wave  function  sufficiently  that  a  plane  wave  solution 
is  inadequate.   The  DWBA  assumes  phase-shifted  spherical  waves  for  the  in- 
cident and  scattered  electron  wave  functions.   Unlike  the  plane  wave  Born 
approximation,  a  form  for  the  charge  distribution  function  must  be  assumed 
before  calculations  are  begun. 

For  elastic  scattering,  the  cross  section  may  be  calculated  using  a 
Fermi  charge  distribution.   However,  this  is  not  appropriate  for  inelastic 
scattering,  in  which  the  outer  nucleons  only  would  be  expected  to  take  part 
in  a  nuclear  transition.   A  good  approximation  was  suggested  by  Tassie 
JRef.  Ill,  in  which  the  transition  charge  distribution  is  assumed  to  be 
the  derivative  of  the  Fermi  function.   Computer  programs  for  calculation 
of  the  differential  cross  sections  ({^)El  and  (^)I1L  have  been  devised  by 
Tuan,  Wright,  and  Onley  \jie£ .    ll]  . 

D.   IDENTIFICATION  OF  EXCITATION  MECHANISMS 

The  form  factor  squared  for  an  inelastic  transition  is  defined  by 

"  (dnL.,E:(dn)m0H'F""u,R- 

Using  the  DWBA,  the  theoretical  variation  of  F^  with  scattering  angle  can 
be  calculated  and  plotted  for  each  EL  and  ML  at  any  incident  energy. 

In  actuality,  one  needs  to  consider  only  L=0  through  about  L=3  since 
higher  terms  are  expected  to  contribute  a  neglible  amount  to  the  inelastic 
scattering.   Each  observed  peak  of  the  inelastic  spectrum  can  then  be 
identified  by  plotting  the  experimental  form  factor  variation  with  angle 
at  each  incident  energy  used  in  the  experiment,  and  comparing  with  the 
theoretical  curves. 
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III.   EXPERIMENTAL  EQUIPMENT  AND  PROCEDURE 

A.  NAVAL  POSTGRADUATE  SCHOOL  LINEAR  ACCELERATOR 

The  NPS  120  MeV  LINAC  was  described  in  detail  by  Barnett  and  Cunneen 
[_Ref.  13J  .   It  is  an  electron  accelerator  similar  in  design  to  the  Stanford 
University  Mark  III  linear  accelerator  described  by  Chodorow  e_t.  al. 
[Ref.  IaJ  .   The  30  foot  NPS  LINAC  is  powered  by  three  klystrons.   Each 
klystron  can  yield  up  to  22  megawatts  of  peak  power.   The  machine  is  pulsed 
at  a  rate  of  60  pulses  per  second  and  each  beam  pulse  is  approximately  one 
microsecond  in  duration.   The  maximum  average  beam  current  is  approximately 
one  microampere. 

Beam  deflection  magnets  and  energy  defining  slits  are  used  to  select 
the  desired  beam  energy,  which  may  range  from  30  to  105  MeV.   The  field  of 
the  beam  deflection  magnet  is  measured  by  a  nuclear  magnetic  resonance 
(NMR)  probe.   Quadrupole  magnets  focus  the  beam  to  a  diameter  of  approxi- 
mately two  millimeters  at  the  target. 

B.  SPECTROMETER  AND  COUNTING  SYSTEM 

The  NPS  LINAC  Spectrometer  was  described  in  detail  by  Oberdier  (_Ref.  15J 
It  is  a  movable,  16-inch  radius,  double-focusing  magnetic  spectrometer 
capable  of  measuring  electrons  scattered  at  angles  from  45   to  150   in 
15  degree  stops. 

The  field  of  the  spectrometer  magnet  is  varied  to  measure  the  energy 
distribution  of  the  scattered  electrons.   Electrons  with  energies  within 
±1.5  per  cent  of  the  spectrometer  energy  setting  are  detected  by  one  of 
ten  NE-102  plastic  scintillation  counters  located  on  the  focal  plane  of 
the  spectrometer  magnet.   Each  detector  serves  as  the  source  of  a  counting 
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channel.   Each  counter  has  a  width  which  allows  detection  in  an  energy  spread 
of  0.3  per  cent  of  the  central  spectrometer  energy.   Two  large  detectors, 
each  of  which  spans  the  same  field  as  the  ten  smaller  detectors,  serve  as 
dual  coincidence  backing  counters.   A  charged  particle  passing  through  a 
front  detector  and  both  backing  detectors  will  register  as  one  count  on  a 
Jorway  Model  1836  Dual  Channel  Scaler  located  in  the  control  room.   The 
counter  gating  circuit  is  synchronized  with  the  accelerator  gun  pulse  so 
that  those  triple  coincidences  that  occur  between  pulses  are  not  counted. 

The  scalers  are  capable  of  recording  periodic  pulses  up  to  150  MHz. 
However,  since  the  beam  is  pulsed,  the  count  rate  was  limited  to  50  counts 
per  second  to  avoid  significant  count' rate  losses.   A  Faraday  cup  is  used 

to  calibrate  the  charge  collected  on  the  Seconday  Emission  Monitor  (SEM) . 

197 
Three  pure    Au  foil  targets  were  used  with  thicknesses  of  two  mils, 

five  mils,  and  ten  mils.   The  targets  were  placed  in  an  evacuated  chamber 
with  1.0  mil  Havar  windows  positioned  to  allow  the  incident  electrons  to 
enter  and  the  scattered  electrons  to  exit.   Since  the  scattering  cross  sec- 
tions decrease  at  the  backward  angles,  five  and  ten-mil  targets  were  used 
in  order  to  increase  the  count  rates  and  thereby  reduce  the  effect  of  sta- 
tistical  fluctuations.   The  two-mil  target  is  recommended  for  future  in- 
vestigations to  maintain  adequate  resolution  (the  ten-mil 

target  gave  an  elastic  peak  width  of  19.5  per  mil  while  the  two-mil  target 
gave  an  elastic  peak  width  of  6.7  per  mil.) 

C.   DATA  ACCUMULATION 

During  each  experimental  "run"  the  beam  energy  and  scattering  angle 
were  held  constant.   The  number  of  counts  in  each  of  the  ten  counting  chan- 
nels and  the  integration  capacitance  and  voltage  were  recorded  automatically 
at  each  spectrometer  energy.   The  spectrometer  energy  was  decreased  in  steps  of 
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0.1  MeV .   The  integration  capacitance  and  voltage  were  chosen  to  maintain  a 
high  count  rate  not  to  exceed  50  counts  per  second.   The  beam  current  was 
decreased  when  measuring  the  elastic  peak  to  avoid  an  excessive  count  rate. 
The  spectrometer  energy  was  reduced  from  100  MeV  to  about  two  MeV  above  the 
elastic  peak.   Elastic  scattering  data  were  taken  at  0.1  MeV  steps  to  about 
four  MeV  below  the  elastic  peak.   When  the  elastic  data  were  accumulated, 
the  beam  was  retuned  to  give  maximum  current  allowed  by  count  rate  restric- 
tions, and  the  spectrometer  energy  was  set  back  several  MeV  to  overlap  the 
elastic  data.   In  order  to  insure  that  the  magnet  always  followed  the  same 
hysteresis  loop,  the  spectrometer  magnet  was  returned  to  a  field  value  cor- 
responding to  100  MeV  and  slowly  decreased  to  the  desired  energy.   It  was 
not  necessary  to  consider  relative  counter  efficiencies  in  this  work.   The 
KIK  code  [Ref.  16]  automatically  corrects  the  energy  of  each  counter  to 
correspond  to  the  true  energy  being  detected. 
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IV.   DATA  ANALYSIS 

A.   DATA  REDUCTION 

The  data  were  collected   on  a  Model  33  teletype  on  paper  and  punched 
tape  simultaneously.   The  punched  tape  was  transcribed  onto  computer  data 
cards  compatible  for  use  on  the  IBM  360-C.   The  organization  of  the  data 
cards  is  described  in  Appendix  A-2.   Fortran  program  DATA  LIST  was  used  to 
print  out  the  transcribed  data  cards  for  verification  against  the  teletype 
"hard  copy."  After  verification  and  correction,  if  required,  the  data 
cards  were  processed  by  fortran  program  KIK.   This  computer  program  assimi- 
lates the  ten  spectra  (one  spectrum  from  each  of  the  ten  spectrometer  chan- 
nels) into  one  total  spectrum,  computes  the  statistical  error,  and  outputs 
the  results  in  both  hard   copy  and  punched  card  form.   The  punched  cards 
from  KIK  are  used  directly  by  fortran  program  SPECTRUM  to  draw  the  total 
spectrum  as  a  point  plot  on  the  CALCOMP  plotters  at  the  NPS  Computer  Center, 
These  spectra  are  displayed  as  Figures  IV-I  through  IV-5. 

The  spectrum  cards  from  KIK  were  then  processed  by  fortran  program  NAW 
[Ref.  1 6 J .   This  program  best  fits  an  analytic  expression  to  the  elastic 
peak  data  and  computes  the  area  under  it  -  the  integration  is  terminated  at 
AE  =  2rs,wherelg  is  the  full-width  of  the  elastic  peak  measured  at  half- 
maximum.   NAW  then  makes  the  appropriate  radiation  corrections  to  the  area, 
as  detailed  in  Section  IV-B.   Experience  has  indicated  that  the  best  value 
for  parameter  L^  is  L,=3  (see  Appendix  A-4). 

Program  NAW  then  determines  an  analytic  expression  to  represent  the 
background  radiation  due  to  the  elastic  tail  in  order  to  make  the  requisite 
subtractions.   The  necessity  for  these  subtractions  will  be  discussed  later 
in  Section  IV-C.   Experience  of  experimenters  at  Institut  Fur  Technische 
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Kernphysik  at  Darmstadt  led  to  the  choice  of  parameter  Lo=3  and  Li^=5 
resulting  in  the  "best"  representation  of  the  radiation  tail.   The  radi- 
ation tail  and  any  constant  background  (.measured  above  the  elastic  peak) 
are  subtracted  from  the  data,  and  program  NAW  then  determines  a  best  fit 
analytic  expression  for  the  inelastic  resonance  peaks.   The  technique  used 
was  to  fix  the  positions  of  the  resonance  peaks  based  on  viewing  the  various 
spectra,  select  peak  widths  (full-width  at  half -maximum)  based  on  other 
published  experimental  results,  and  then  dynamically  adjust  these  parameters 
until  one  best  fit  for  all  experimental  spectra  was  obtained.   Program  NAW 
displays  these  results  both  graphically  and  in  tabular  form.   Furthermore, 
it  integrates  the  analytic  expression  for  the  inelastic  resonances  to  deter- 
mine the  areas  under  the  peaks  and  displays  the  ratio  of  the  area  of  the 
inelastic  peak  to  the  area  of  the  elastic  peak  for  each  resonance. 

Using  the  incident  electron  energies  calculated  by  NAW,  fortran  program 
MARPHASE  uses  the  Rawitscher-Fischer  Partial  Wave  Code  [Ref.  8  J  to  calcu- 
late the  elastic  scattering  parameters;  the  results  are  displayed  in 
Table  IV- 1.   Since  the  differential  scattering  cross  sections  are  propor- 
tional to  the  area  under  the  resonance  peaks  one  can  easily  determine  the 

inelastic  scattering  cross  sections,  or  more  usefully,  the  square  of  the 

o 
inelastic  scattering  form  factor,  F  : 

p2  _  (c/o-/dJi)inei  _  Ai  >  (d<r/d.si)elas 
(d<r/dJlW  "  Ae  '  (d<r/dn)mott  . 

The  inelastic  form  factors  squared  may  now  be  graphed  as  a  function  of  the 
scattering  angle  and  parametized  by  incident  energy.   To  identify  the  ex- 
citation involved,  one  may  compare  these  angular  plots  with  the  theoretical 

values  as  determined  by  a  DWBA  code.   Fortran  program  DUKE  [Ref.  12j  uses 

2 
a  DWBA  code  to  calculate  F   for  both  electric  and  magnetic  transitions. 
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TABLE  IV- 1 


ELASTIC  SCATTERING  PARAMETERS  FOR 
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AU 


0 

E{(MeV) 

da/da 

(do /da) 

F2 

75° 

75.15 

2.52xl0-1 

2.62 

9.59xl0~2 

75° 

90.19 

7.83xl0"2 

1.82 

4.30x10-2 

90° 

75.20 

5.62xlO-2 

1.14 

4.91x10-2 

90° 

90.19 

2.39xl0_2 

7.95x10   l 

3.00x10-2 

120° 

75.06 

8.75xl0-3 

2.55x10-! 

3.43x10-2 

Differential  cross  sections  are  in  units  of  fib/sr. 

A  transition  strength,  B,  of  unity  is  assumed  for  all  transitions  in  these 

calculations.   While  the  output  labeled  as  F0RM  is  the  theoretical  value 

o  2 

for  F^  for  electric  transitions,  the  corresponding  output  is  F  multiplied 

by  the  factor  V-r/V^  (labeled  QVLT  on  the  output)  for  magnetic  transitions. 

Once  the  theoretical  values  for  F  have  been  obtained  the  experimental 

values  for  the  transition  strength  can  be  calculated  by  simply  averaging 

the  values,  i.e. 

b4S(h:/i-). 


Fortran  program  BCALC  computes  B  and  then  multiplies  the  theoretical  values 

2  2 

of  F  by  B.   These  are  graphed,  and  the  experimental  values  of  F  are  then 

plotted  on  the  same  graph.   This  procedure  is  done  with  each  identified 
resonance  for  all  reasonably  possible  transitions.   The  transition  type 
identification  is  made  by  determining  which  graph  exhibits  the  best  fit  of 
the  experimental  points  to  the  theoretical  curves  and  also  yields  a  reason- 
able value  for  B.   To  avoid  subjective,  and  thereby  possibly  prejudiced, 
evaluations  of  the  best  fit,  program  BCALC  computes  the  normalized  quanti- 
ties, VAR,  which  are  compared  to  determine  the  best  fit  -  the  smaller  the 


value  of  VAR,  the  better  the  fit, 
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VAR=Z(^£H£). 


These  values  may  be  found  in  Table  V-2. 


B-Ft- 


B.   RADIATION  CORRECTION  CALCULATIONS 

There  are  processes  occurring  during  elastic  scattering  that  result  in 
low  energy  electrons  which  constitute  the  so-called  elastic  radiation  tail. 
Since  it  is  not  practicable  to  integrate  the  electron  spectrum  under  the 
elastic  peak  to  zero  energy,  it  is  conventional  to  apply  three  correction 
factors.   It  is  sufficient  to  correct  for  such  processes  as  follows: 


Bsb  Stw«  • 


where  Kg,  Kg,  and  K,  are,  respectively,  the  Schwinger,   Bremsstrahlung,  and 
Electron-Electron  Scattering  corrections. 
1 .   Schwinger  Correction 

The  Schwinger  term  is  by  far  the  most  significant  of  the  three 
corrections  being  discussed.   This  term  corrects  for  virtual  photon  emission 
and  reabsorption  and  real  photon  emission  occurring  in  the  field  of  the  tar- 
get nucleus  -  both  of  which  result  in  large-angle  scattering  of  the  incom- 
ing electron,  i.e.,  during  the  scattering  event.   Schwinger's  original 
rk  [Refs.  17,  18]  gives  the  relation 

';    J-L    ;     If  ({))    sin  'Of-) 


WO 


Ks  ■   exp[-^/rin-f3--Hl.hnflIa.sin  9/2\      _1 

I  n  VL    ae     12J  L   W*  '    2 
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where  a  is  the  Sommerfield  fine  structure  constant,  E»  is  the  most  probable 
energy  of  the  outcoming  electron,  Ae  is  the  width  of  the  integration  (in 
program  NAW  AE  is  set  at  twice  the  elastic  peak  full-width  at  half-maximum), 
and  f(0)  is  the  Spence  Function.   Many  investigations  have  since  improved 
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upon  this  formulation,  notably  the  work  of  Tsai  [_Ref.  1 9  J  and  that  of 
Meister  and  Yennie  [^Ref.  20J  . 

2.   Bremsstrahlung  Correction 

This  term  corrects  for  the  energy  radiated  via  bremsstrahlung  during 
small-angle  deflections  caused  by  the  electron's  interactions  with  the  fields 
of  nuclei  other  than  the  elastic  scattering  center  both  before  and  after  the 
elastic  scattering  event.   This  energy  loss  mechanism  was  first  studied  by 
Bethe  and  Heitler  [Ref.  2lj  ,  and  the  correction  factor  for  this  mechanic 
was  given  by  Hofstadter  |Ref.  22J  as 


Lsm 


K»=eXp["  X^2-|nlt]   ' 


where  X0  is  the  radiation  length  and  t  is  the  target  thickness,  both  in 
g/cm  .   Program  NAW  uses  the  formalism  of  Barber,  et_.  aT.   Ref.  23   to  give 


K. 


It/     1      E. 
P[2cos  9/2\X0(Ei-Ef)  E. 


( V  +  Ef  2>  -  g   + 

E.E 


'if 


3     91n(183exp{-lnzv 


*»>])!. 


3.   Electron-Electron  Scattering 

This  term,  sometimes  called  the  Moller  term,  is  almost  negligible. 
The  term  corrects  for  electron-electron  scattering  in  the  target  that  occurs 
before  or  after  the  elastic  scattering  event.   The  effect  is  to  produce  low- 
er energy  electrons,  almost  in  the  original  direction,  similar  to  the  effect 
of  the  small-angle  bremsstrahlung.   By  Barber's  formalism  we  have 


KM  =  exp< 


2nz 


] 


2E_L  +  1 


1 


1 


Sf"   Ef(E.-Ef)(E.+l)2  +  (E±-Ef)^  +  (E.+l)^ 


C.   RADIATION  TAIL  CALCULATIONS 

When  calculating  the  radiation  corrections  we  were  concerned  only  with 
determining  how  many  elastically  scattered  electrons  were  degraded  to  an 
energy  beyond  the  integration  cut-off,  Ae.   The  calculation  of  the  radiation 
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tail  is  concerned  with  the  energy  distribution  of  electrons  lost  from  the 
elastic  peak.   Program  NAW  has  the  option  of  using  either  Schiff's  method 
[Ref.  24J  or  the  more  recent  method  developed  by  Ginsberg  and  Pratt  [Ref.  25j. 
The  latter  was  used  in  this  experiment  resulting  in  an  approximation  to  the 
true  radiation  tail  -  while  the  approximation  is  crude,  it  does  have  the 
proper  slope.   NAW  then  best  fits  a  fourth-degree  polynomial  to  the  cal- 
culated values  (output  in  program  NAW  as  SSW) .   The  resulting  expression 
for  SSW  is  then  fit  to  the  raw  data  spectrum  after  any  constant  background 
counts  have  been  subtracted.   Of  the  11  fit  options  available,  the  one  used 
in  this  work  was 

UNT  =  A-SSW  +  B  +  OX, 
where:   UNT  is  the  portion  of  the  measured  counts//!  C  attributed  to  the 
elastic  radiation  tail;  A,  B,  and  C  are  numerical  fitting  parameters  deter- 
mined by  NAW;  and  X  =  E-E-g  (Eg  is  a  constant  determined  by  NAW)  .   The  radi- 
ation tail  obtained  is  then  subtracted  from  the  inelastic  data  to  yield 
the  unfolded  spectrum. 

D.   ERROR  ANALYSIS 

1 .   Systematic  Errors 

In  this  experiment,  the  ratio  of  the  inelastic  to  the  elastic  cross 

section  is  computed.   Since  the  inelastic  and  the  elastic  data  were  always 
taken  during  the  same  experimental  run,  most  systematic  errors  cancelled. 
These  would  include  errors  due  to  variations  in  solid  angle,  scattering 
angle,  SEM  efficiency,  target  thickness,  target  purity,  and  absolute  spec- 
trometer efficiency.   On  several  occasions  after  the  inelastic  data  points 
were  collected,  the  elastic  peak  points  were  repeated  to  verify  that  the 
number  of  counts  per  microcoulomb  remained  constant,  within  statistical 
variations,  during  a  run. 
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2.   Random  Errors 

The  statistical  fluctuations  of  the.  number  of  electrons  scattered 
for  a  given  number  of  incident  electrons  obeys  Poisson  statistics,  in  which 
the  standard  deviation  is  given  by  SN(E)  =  J  N(E) ' ,  where  N(E)  is  the  total 
number  of  counts  at  scattered  energy  E.   The  fractional  error  is  therefore 
given  by: 

\/N(E)' 

Fractional  error  =  . 

N(E) 

Program  KIK  computes  the  energy,  E,  for  each  counter  at  each  spectrometer 
setting,  sums  the  ten  spectra  to  yield  one  total  spectrum,  and  tabulates 
the  total  counts  per  microcoulomb  and  the  associated  fractional  error. 
This  error  varied  from  one  percent  to  three  percent.   The  effect  of  this 
error  on  the  final  form  factors  is  calculated  by  NAW  as  follows.   Let  the 
total  counts  detected  in  an  energy  region  encompassing  the  width  of  an  in- 
elastic peak  be  N™  and  the  counts  attributed  to  the  inelastic  peak  be  N, 
(thus  N  -  N   is  attibuted  to  radiation  tail  and  background  in  this  region) , 
Then  the  fractional  error  in  the  cross  section  is  given  by   JN^' /Nj.  Due 
to  the  varying  strength  of  each  transition  and  the  fact  that  the  total  num- 
ber of  counts  was  lower  for  certain  runs,  these  errors  varied  greatly  in 
size.   They  are  tabulated  in  Table  V-l. 

Other  random  errors  include  changes  in  the  beam  current,  variation  in 
SEM  efficiency  during  a  run,  energy  shifts  of  the  incident  beam,  leakage 
current  in  the  beam  current  integrators,  error  in  measurement  of  the  back- 
ground, uncertainties  in  the  position  of  the  peaks,  and  uncertainties  in 
the  widths  of  the  peaks. 

Variations  in  the  beam  current  affect  the  values  obtained  for  the 
cross  sections  since  count  rate  losses  increase  with  current.   Although 
these  losses  were  initially  believed  small  for  count  rates  up  to  50  counts 
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per  second,  current  studies  indicate  that  this  is  probably  not  so.   The 
magnitude  of  this  error  is,  as  yet,  unknown,  but  may  well  contribute  as 
much  as  an  additional  25  percent  error  to  the  experimental  values  of  the 
form  factors. 

The  possibility  of  variations  of  the  SEM  efficiency  during  a  run  was 
not  investigated  in  this  experiment.   Measurements  by  Berard  and  Traverso 
|Ref .  26j  show  that  the  SEM  efficiency  remains  constant  within  0.09  percent. 

Variations  in  line  voltage  cause  fluctuations,  not  only  in  beam  energy, 
but  in  spectrometer  and  beam  deflection  magnet  fields  as  well.   The  effect 
is  to  shift  some  electrons  from  one  counter  to  another.   However,  the  elec- 
tronic regulation  of  the  spectrometer  and  deflection  magnets  reduce  this 
error  to  a  negligible  level.   Leakage  currents  in  the  beam  integration  cir- 
cuit were  measured  by  Topping  JRef.  27  J  and  were  found  to  be  negligible 
also. 

The  background  contributes  less  than  1  percent  of  the  area  of  the  elas- 
tic peak,  so  that  the  fractional  error  due  to  the  background  subtraction 
is  negligible  -  particularly  when  compared  to  the  uncertainty  involved  in 
the  subtraction  of  the  radiation  tail. 

The  error  involved  in  fitting  inelastic  peaks  to  the  data  is  difficult 
to  assess.   The  experimental  form  factors  calculated  by  program  NAW  depend 
on  the  choice  of  location  and  width  of  each  of  the  inelastic  peaks.   It  is 
estimated  that  these  uncertainties  contribute  a  20  percent  error  to  the 
experimental  values  of  the  form  factors. 

The  total  fractional  error  for  any  given  calculation  of  an  experimental 
form  factor  can  now  be  obtained  as  follows: 


s=/siF, 


allMwrcts 
•♦  trrt  r 
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where  S-  is  the  fractional  error  due  to  the  i'th  mechanism.   Since  three 

factors  dominate  the  others,  the  fractional  error  can  be  approximated  by 

I     2        2     2 
S  =   .  (.25)   +  (.2)   +  S   ,  where  S   is  the  fractional  error  due  to  Poisson 

counting  statistics  (found  in  Table  V-l).   In  most  cases,  the  total  esti- 
mated error  was  approximately  40  percent. 
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V.   EXPERIMENTAL  RESULTS  AND  CONCLUSIONS 

A.   EXPERIMENTAL  RESULTS 

The  results  of  five  experimental  runs  are  reported  here.   Three  other 
runs  were  conducted  but  are  not  reported.   Two  runs  conducted  at  a  scatter- 
ing angle  of  150  had  counts  that  were  so  small  that  any  resonances  present 
are  masked  by  the  statistical  fluctuations.   One  run  conducted  at  an  inci- 
dent energy  of  90  MeV  and  a  scattering  angle  of  120  using  a  ten-mil  target 
had  insufficient  resolutions  to  extract  any  structure  from  the  resultant 
spectrum.   Prior  to  presenting  the  experimental  results,  a  brief  description 
of  the  procedure  by  which  program  NAW  was  used  follows. 

The  initial  positions  of  the  resonances  were  estimated  by  viewing  the 
experimental  spectra,  resulting  in  the  values  of  7.4,  8.8,  10.6,  13.8  and 
14.2,  and  23.0  MeV.   Since  the  fit  by  NAW  for  the  radiation  tail  is  only 
considered  valid  for  a  range  of  about  20  MeV  and  the  theoretical  values  for 
the  radiation  tail  are  known  to  become  progressively  inaccurate  with  in- 
creasing excitation  energy,  the  23.0  MeV  resonance  was  not  investigated  in 
detail  (preliminary  results  indicate  that  this  resonance  is  probably  an 
E2  transition).   The  remaining  four  resonances  were  given  to  program  NAW, 
with  widths  of  3.0,  2.5,  3.0,  and  4.75  MeV  respectively,  and  the  program 
was  allowed  to  vary  the  positions  of  the  last  two  resonances.   Each  run 
resulted  in  somewhat  different  values  for  the  positions  of  the  resonances, 
and  average  positions  were  determined.   The  averaged  values  were  then  en- 
tered into  program  NAW  and,  keeping  all  four  resonances  fixed  in  energy, 
the  widths  were  all  allowed  to  be  determined  by  NAW.   At  this  stage,  for 
each  resonance,  the  average  value  of  the  width  for  all  the  runs  was  calcu- 
lated.  Entering  these  averaged  values  for  the  widths  into  NAW,  the  7.4 
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and  8.8  MeV  positions  were  allowed  to  vary  and  again  the  averaging  process 
was  followed.  The  new  positions  were  given  to  NAW,  the  widths  were  once 
more  allowed  to  vary,  and  the  results  were  again  averaged.  The  final  re- 
sults of  this  procedure  were  resonances  at  7.14,  9.12,  10.9,  and  13.8/14.2 
MeV  with  respective  widths  of  3.13,  2.89,  3.06,  and  4.89  MeV.  The  peak  in 
the  vicinity  of  14  MeV  was  found  to  be  incident  energy  dependent.  At  an 
incident  energy  of  75  MeV  the  peak  maximum  was  located  at  13.8  MeV,  and  at 
90  MeV  the  peak  maximum  was  located  at  14.2  MeV.  Hereafter  this  resonance 
will  be  referred  to  as  the  14  MeV  peak. 

Program  NAW  determined  an  analytic  expression  to  best  fit  the  data, 
resulting  in  a  smoothed  spectrum.   In  also  determines  an  analytic  expres- 
sion for  the  radiation  tail;  these  functions  are  displayed  graphically  in 
Figures  V-l  through  V-5.   The  radiation  tail  and  any  constant  background 
were  then  subtracted  from  the  smoothed  spectrum,  and  analytic  expressions 
for  the  four  inelastic  resonances  were  determined  assuming  a  Breit-Wigner 
form.   The  areas  under  the  inelastic  resonances  were  then  computed  and  the 
ratio  of  each  of  these  areas  to  the  area  under  the  elastic  peak  was  calcu- 
lated along  with  the  fractional  error  due  to  counting  statistics  (Poisson) . 
Multiplying  A-/A  by  the  square  of  the  elastic  scattering  form  factor  (cal- 
culated by  program  MARPHASE  and  displayed  in  Table  IV-1)  gives  the  square 
of  the  inelastic  form  factor.   The  results  of  these  calculations  are  dis- 
played in  Table  V-l. 

2 

Using  the  experimental  values  of  F  at  the  different  angles  and  ener- 
gies and  the  theoretical  values  as  calculated  by  DUKE  as  input  data,  program 
BCALC  determines  the  transition  strengths,  B,  and  the  values  of  VAR  for 
the  various  transitions.   It  also  prints  out  (for  subsequent  plotting) 
experimental  and  theoretical  values  of  F  based  on  the  calculated  value 
of  B.   The  numerical  results  are  displayed  in  Table  V-2  and  the  graphs  are 

reproduced  as  Figures  V-6  through  V-9. 
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TABLE  V-l 


EXPERIMENTAL  INELASTIC  SCATTERING  PARAMETERS  FOR    AU 


197 


Ex  (MeV) 

t(mi|) 

0 

EQ(MeV) 

Aj/Ae 

Percent 
Error 

F2 

2 

75° 

75.04 

2.24xl0-2 

24 

2.14xl0"3 

2 

90° 

75.07 

2.41xl0~2 

30 

1.18xl0"3 

7.14 

5 

120° 

74.82 

3.45x10-2 
3.54x10  l 

12 

1.18xl0-3 

2 

75° 

90.07 

17 

1.52xl0~3 

2 

90° 

90.05 

2.03xlO"2 

24 

6.10xl0"4 

2 

75° 

75.04 

7.07xl0-3 

18 

6.78x10-4 

2 

90° 

75.07 

7.05xl0-3 

24 

3.46xl0-4 

9.12 

5 

120° 

74.82 

4.40xl0-3 

41 

1.51xl0-4 

2 

75° 

90.07 

1.20xl0-2 

11 

5.16xl0-4 

2 

90° 

90.05 

6.93xl0~3 

22 

2.08xl0"4 

2 

75° 

75.04 

1.06xl0"2 

9 

1.02xl0"3 

2 

90° 

75.07 

l.lOxlO-2 

13 

5.38xl0-4 

10.90 

5 

120° 

74.82 

2.08xl0~3 

64 

7.14xl0-5 

2 

75° 

90.07 

1.37xl0-2 

9 

5.90xl0-4 

2 

90° 

90.05 

8.79xl0-3 

16 

2.64xl0~4 

2 

75° 

75.04 

1.56xl0-2 

7 

1.50xl0-3 

? 

90° 

75.07 

1.32xl0-2 

14 

6.49xl0"4 

14.00 

5 

120° 

74.83 

5.69xl0-4 

298 

1.95xl0-5 

2 

75° 

90.07 

1.59x10-2 

10 

6.84xl0~4 

2 

90° 

90.05 

l.llxlO-2 

15 

3.32xl0~4 

TABLE  V-2 
MULTIPOLARITY  PARAMETERS 


Ex(MeV) 

Type 

B(e*fmlL) 

VAR 

7.14 

Ml 

163.9 

0.956 

7.14 

M2 

14,895.9 

0.511 

7.14 

E2 

31,655.0 

0.880 

9.12 

Ml 

48.4 

1.574 

9.12 

E2 

7,819.7 

0.125 

10.90 

El 

91.2 

1.105 

10.90 

E2 

8,822.9 

0.582 

14.00 

El 

89.9 

1.502 

14.00 

E2 

10,527.6 

1.390 
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B.   CONCLUSIONS 

Transition  assignments  were  made  as  discussed  in  Section  IV-A.   The 
7.14  MeV  resonance  appears  to  be  an  M2  transition,  but  this  assignment  can- 
not be  made  conclusively  due  to  lack  of  scattering  data  at  the  background 
angles.   Future  investigations  should  concentrate  on  obtaining  data  at  the 
backscattering  angles,  since  the  cross  sections  for  magnetic  transitions 
increase  at  these  angles.   It  may  also  be  useful  to  obtain  data  at  the  small 
forward  angles,  in  order  to  move  away  from  the  peaks  in  the  form  factors. 
The  heretofore  unreported  resonance  at  9.12  MeV  is  concluded  to  be  an  E2 
transition.   The  assignment  of  the  10.9  MeV  resonance  as  an  E2  is  in  agree- 
ment with  previous  unpublished  work  at  Darmstadt  fRef.  4_| .   The  results  of 
our  work  indicate  that  the  resonance  collectively  referred  to  as  14  MeV 
shifts  with  incident  energy.   This  shift  has  been  explained  by  Pitthan 
[Ref .  2j  ,  who  assumed  that  the.  broad  resonance  in  the  cross  section  is  cal- 
culated from  a  broad  nuclear  excitation  strength  function.   This  strength 
function  varies  primarily  with  excitation  energy,  but  also  varies  with  in- 
cident beam  energy  and  scattering  angle.   Thus  the  maximum  of  the  inelastic 
peak  may  shift  when  beam  energy  is  changed.   The  results  of  this  work, 
coupled  with  the  results  of  the  Darmstadt  group  f"  4J  (conducted  at  lower 
incident  energies)  indicate  this  resonance  is  an  El  giant  resonance. 

The  results  described  above  are  recapitulated  in  Table  V-3.   All  pre- 
viously reported  Ref.4|  resonances  appeared  in  our  spectra.   The  observed 
strengths  \i7ere  not  in  exact  agreement  with  the  previous  results,  but  dis- 
agreement can  be  explained  by  observed  deficiencies  in  the  counting  system. 
Furthermore,  two  additional  resonances  were  observed:   9.1  MeV,  with 
B(E2)  =  S.OOOe2^;  and  23  MeV  which  is  probably  an  E2  transition.   The 
limitations  of  the  determination  of  the  elastic  radition  tail,  and  insuf- 
ficient experimental  data  precluded  detailed  investigations  of  the  latter 
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resonance.   We  have  demonstrated  that  giant  resonances  can  be  measured  at 
the  NPS  LINAC  Facility.   The  various  experimental  techniques  developed,  the 
problem  areas  identified  and  rectified  where  possible,  and  the  computer 
data  handling  procedures  developed  will  enable  future  investigations  of 
giant  resonances  by  the  LINAC  group  to  proceed  much  more  rapidly  and  accu- 
rately.  Future  investigators  should  spend  a  fair  amount  of  time  developing 

new  techniques  for  fitting  the  elastic  radition  tail,  because  its  effect 

2 

on  the  resulting  values  for  F   is  paramount. 


TABLE  V-3 
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IDENTIFIED  GIANT  MULTIPOLE  RESONANCES  IN    AU 


Ex  (MeV) 

Q  (MeV) 

Traruition 
Type 

B(ezfm*L) 

7.1 

9.1 

10.9 

14.0 

3.1 
2.9 
3.1 
4.9 

M2 
E2 
E2 
El 

15,000 

8,000 

9,000 

90 
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APPENDIX  A 
COMPUTER  PROGRAMS 

Seven  FORTRAN  programs  were  used  to  perform  the  necessary  calculations 
for  this  work.   The  computer  programs  are  reproduced  on  the  following  pages. 
Each  program  listing  will  be  prefaced  by  a  brief  explanation  of  what  it  does, 
and  a  parameter  list  with  explanatory  remarks  (program  options  not  directly 
applicable  to  this  and  related  works  will  not  be  listed).   Additionally,  if 
any  special  computer  control  cards  are  required  to  execute  the  program  at 
the  NPS  Computer  Facility,  a  sample  FORTRAN  deck  set-up  will  be  included. 

A.  DATA  LIST 

DATA  LIST  prints  the  raw  data  cards  in  the  teletypewriter  format  used 
by  the  NPS  LINAC  Data  Acquisition  System,  segregating  by  integral  data  sets, 
so  as  to  aid  in  verifying  the  accuracy  maintained  during  the  processing  from 
punched-tape  output  to  punched-card  output.   This  program  requires  only  one 
parameter  card: 

TEST,  TEXT  (A4,19A4) 

TEST      The  first  four  characters  on  the  first  data  card. 
Usually  a  blank  followed  by  three  digits. 

TEXT      Any  alphanumeric  message  to  be  printed  on  the  first 
page  of  the  data  listing. 

B.  KIK 

KIK  can  process  up  to  499  data  sets,  each  set  consisting  of  the  scaler 
readings  for  each  of  the  ten  spectrometer  channels  of  the  NPS  LINAC  Count- 
ing System.  KIK  compiles  the  ten  separate  spectra  into  one  total  spectrum 
and  calculates  the  errors,  giving  the  results  on  both  hard  copy  and  punched 
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cards.   A  FORTRAN  Deck  set-up  consists  of  the  requisite  control  cards,  four 
parameter  cards,  the  data  deck,  and  finally  two  trail  cards. 

1.  TEXT       Any  identifying  alphanumeric  remarks  desired    (20A4) 

2.  K,  FRE0,  FREU,  DF ,  TAU,  FEA,  RH0,  THETA  (I10,7F10. 5; 

K=l        Program  control  parameter  signifying  the  start 
of  a  spectrum. 

FRE0       Upper  spectrometer  energy  limit  for  the  spectrum 
in  MeV. 

FREU       Lower  spectrometer  energy  limit  for  the  spectrum 
in  MeV. 

DF         Spectrometer  energy  step  width  in  MeV.   KIK  will  in- 
terpolate between  missing  data  sets. 

TAU*       Counting  system  dead  time/pulse  length. 

FEA*       Incident  accelerator  energy  (MeV) . 

RH0*       Target  thickness  (g/cm2). 

THETA*     Electron  scattering  angle  (degrees) . 
All  parameters  followed  by  an  asterisk  are  for  identification  purposes  only 
and  are  not  used  in  the  calculations.   They  may  be  left  blank  if  desired. 

3.  a(l)...a(10)  (1014) 

Spectrometer  counting  system  normalization  factors 
relative  to  counter  number  5.   The  a.'s  are  defined  by: 


a .  = 

x 


ENERGY (5) 
ENERGY (i) 


*10' 


Input/Output  Control 

LL,  LU,  L0,  LF0R,  INV,  FKXAD,  FKEN 

LL=1 


(512.2F10.5) 


LU,  L0     The  spectrum  of  counter  I  specified  by  LU<KL0  will  be 
punched  out  onto  cards. 

LF0R=O 

INV>0      The  order  of  the  input  data  is  reversed;  i.e.,  from 
lower  to  higher  energies. 
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FKXADX)        Tne  input  capacitor  voltage  (full-scale  setting) 
will  be  scaled  to  units  of  lO-^  volts  by  internal 
multiplication  by  this  factor.   For  the  NPS  LINAC 
Facility,  FKXAD=100  with  the  voltage  being  record- 
ed in  millivolts. 

FKEN>0         The  input  energy  will  be  scaled  to  units  of  KeV 
by  internal  multiplication  by  this  factor.   For 
the  NPS  LINAC  Facility,  FKEN=100  with  energies 
recorded  in  MeV. 

The  data  cards  now  follow  with  two  cards  per  data  set.   The  actual  format 

used  is  quite  complex,  since  the  program  reads  data  around  the  colons  and 

reconstructs  decimal  numbers.   The  general  order  of  the  input  data  fields 

is  as  follows: 

Card  1:   E0,  VFS,  E,  C,  T,  V,  DUMMY,  C0(1) . . .C0(3) 

Card  2:   C0(4) . . .C0(1O) ,  B,  A 

E0      Incident  accelerator  energy 

VFS     Full-scale  voltage  setting  on  t-he  capacitor  (mv), 

E       Spectrometer  energy  setting  (MeV) . 

C       Canacitance  of  the  collecting  caDacitor  QaF) - 

T       Elapsed  time  (sec) . 

V       Final  voltage  on  the  capacitor  (mv) . 

DUMMY  =  000000 

Spectrometer  counter  readings. 

Spectrometer  backing  counter  reading. 

Spectrometer  'accidentals'  counter  reading. 


C0CD 

B 
A 
5.   Trail  Cards. 


When  a  card  is  encountered  with  columns  1-7  left  blank,  that  card   and 
the  immediately  following  card  are  interpreted  as  control  cards,  with  the 
following  meanings: 

VFS<0  RMIN=VFS*10-3 


=0 


RMIN=10" 
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Two  points  differing  from  each  other  in  energy  by 
less  than  RMIN  will  be  included  as  one  point  on 
the  spectrum. 

C0(I)>O  Counter  I  will  not  be  considered  in  compiling  the 

total  spectrum. 
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C.   SPECTRUM 

SPECTRUM  uses  the  output  cards  from  KIK  as  its  input  data,  displaying 
the  spectrum  as  a  'point  plot'  drawn  on  the  NPS  Computer  Facility's  CALCOMP 
Plotter.   The  graph  will  be  15  inches  long  and  9  inches  high.   SPECTRUM  re- 
quires three  parameter  cards;  the  first  two  cards  contain  alphanumeric  iden- 
tifying text  in  columns  1-48  and  the  third  is  as  follows: 
NPTS,  XSCL,  ESCL,  E0FFST,  X0FFST  (13 ,7X,4F10. 0) 

NPTS      Numer  of  points  to  be  plotted,  with  the  restriction  that 
NPTS>30n  +1;  n  is  any  positive  integer. 

XSCL      Ordinate  scale  in  units/inch  (only  one  significant  digit.) 

ESCL      Abscissa  scale  in  similar  fashion. 

E0FFST    Position  of  the  abscissa  zero  point.   Usually  for  the 

elastic  resonance  E0FFST=Ee-^  +  2;  for  the  inelastic  spectrum 
E0FFST=Ee2  which  results  in  the  points  being  displayed  at 
excitation  energies. 

X0FFST    Position  of  the  ordinate  zero  point. 
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D.   NAW 

NAW  uses  the  output  from  KIK  as  its  input  data.   This  program  calculates 
the  ratio  of  the  area  under  the  inelastic  resonance  peaks  to  the  area  under 
the  elastic  resonance  peak.   Because  of  its  extreme  length  (3134  FORTRAN 
cards)  the  program  listing  will  not  be  reproduced;  documentation  may  be 
found  in  reference  15.   NAW  accepts  up  to  198  elastic  data  points,  498  in- 
elastic data  points,  and  25  inelastic  resonances  per  spectrum. 

1.  Text   Alphanumeric  identifying  remarks.  (20A4) 

2.  EICH,  PN2,  PN3,  PN4,  PN5,  PN6 ,  PN?  (7F10.5) 

EICH   Spectrometer  energy  normalization  factor  to  scale  all 
input  energies  to  MeV. 

PN .    Starting  values  for  the  various  parameters  being  used  to 

fit  the  elastic  spectrum  with  an  analytic  expression.   They 
are  used  only  when  Lj>0  on  card  #7. 

PN2   Half-width  of  a  gaussian  to  fit  the  high  energy  side  of 
the  elastic  peak. 

PNo   Half-width  of  a  gaussian  to  fit  the  low  energy  side  of  the 
elastic  peak. 

PN^   Matching  point  for  the  low  energy  gaussian  and  a  hyperbolic 
segment.   PN,  is  specified  in  MeV  from  the  position  of  the 
zero  point  and  furthermore  is  a  negative  number. 

..  PN,-    Position  of  the  zero  point.   When  PN^=0  the  zero  point  is 
taken  to  be  the  position  of  the  elastic  resonance  peak. 

PNg>0  Height   of  the  second  elastic  resonance  peak. 

PNy>0  Position  of  the  second  elastic  resonance  peak. 

3.  EE,  THETA,  D,  UNT,  GINT  (5F10.50) 
EE    Elastic  resonance  energy  in  MeV. 

THETA  Electron  scattering  angle  in  degrees. 

D     Target  thickness  in  g/cm-. 

UNT   Constant  background  at  the  high  energy  end  of  the  elastic 
resonance  spectrum  in  counts//iC. 
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GINT  Lower  integration  limit  for  calculating  the  area  under  the 

elastic  peak.   This  value  is  only  used  when  L.=0  on  card  #7. 
At  all  other  times  the  lower  limit  is  fixed  at  E  ,  minus  twice 
the  peak  width. 

4.  EEL  ,  EEL  ,  GSSW  ,  GSSW  ,  DSSW,  DWAL  (6F10.5) 

o      u '      O       u '       ' 

EELQ   High  energy  limit  for  fitting  the  elastic  peak. 

EEL   Low  energy  limit  for  fitting  the  elastic  peak. 

GSSWQ  Upper  limit  for  calculating  the  radiation  tail. 

GSSWU  Lower  limit  for  calculating  the  radiation  tail. 

DSSW  Energy  step  to  be  used  in  calculating  the  radiation  tail. 

DWAL   See  L2  on  card  #7. 

5.  E0,  Eu  (2F10.5) 
E     High  energy  limit  for  fitting  the  inelastic  spectrum. 

Eu    Low  energy  limit  for  fitting  the  inelastic  spectrum. 

6.  GRj,  GRl,  GR2,  GR2,  GR3,  GR3  (6F10.5) 

Regions  where  it  is  expected  that  only  background  contributes  to 
the  spectrum;  i.e.,  between  resonances  in  the  inelastic  spectrum. 
There  must  be  at  least  one  set  GRQ,  GRU  input,  with  a  maximum  of 
three  such  sets. 

7.  Ll5  L2,  ...,  L22  (2213) 
Li     Analytic  expression  used  to  fit  the  elastic  resonance.  • 

=0   Simple  curve  sections  (ARC3) . 

=1   A  continuous  curve  consisting  of  txro  gaussians  and  a  seg- 
ment of  a  hyperbola  (ARC5). 

=2   The  elastic  resonance  will  be  fit  as  if  L^=l  and  the  in- 
elastic resonances  will  be  fit  using  the  Breit-Wigner  formula. 
Each  curve  will  be  multiplied  by  a  straight  line  with  equa- 
tion Y=l  +  ¥±(Eyr-Eox)'      pi  will  be  specified  on  card  #10; 
Ex  is  the  excitation  energy,  and  EQX  is  the  position  of  the 
maximum. 

=3   The  elastic  resonance  will  be  fit  as  if  Li=2  except  that 
NAW  will  determine  the  values  P-^  and  card  #10  will  be  left 
blank. 

Lo    Analytic  expression  used  to  represent  the  background  due 
to  the  radiation  tail. 
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=1   BG  =  Pj^SSW 


=2   BG  =  Pi*SSW  +  ?2 

=3   BG  =  Pj*SSW  +  P2   +  P3X 


=4  BG  =  P 
=5  BG  =  P 
=6  BG  =  P 
=7  BG  =  P 
=8  BG  =  P 
=9  BG  =  P 
=10  BG  =  P 
=11  BG  =  P 


+  ssw 


+  p2x  +  ssw 

+  p2x  +  p3x2  +  ssw 


+  P2A  +  SSW 

+  P2A  +  P3a2  +  SSW 

+  P2B  +  P3B2  +  P^X  +  SSW 


+  SSW*(P2  +  P3X) 

+  SSW*(P2  +  P3X  +  P4X2) 

P.  will  be  determined  by  NAW. 
l  J 

SSW  =  Theoretically  calculated  radiation  tail, 


X  =  E  -  E_ 
lated  by  NAW 


B.   The  reference  energy  Eg  will  be  internally  calcu- 


A  =  1-(E  ,  -  E). 
el 

B  =  l/(Eel  -  E  -  DWAL) .  Should  L2=9  be  selected,  then  DWAL  must 
be  input  on  card  #4.  If  L2/9  then  the  space  for  DWAL  on  card  #4 
should  be  left  blank. 

Lo    The  number  of  resonances  with  specified  energies.   If  L-XO 
then  the  resonances  will  be  kept  fixed  relative  to  E  (1), 
but  Ex(l)  will  be  positioned  by  NAW. 

L,     The  number  of  resonances  with  energies  to  be  determined  by 
NAW.   Note  that  L3  +  L^  must  not  exceed  25. 

Lr     The  number  of  line  widths  to  be  determined  by  NAW. 

=0  All  line  widths  will  be  specified  on  card  //9. 

>0   E.g.  =3,  means  that  three  line  widths  are  free  and  in  fact 
are  Ex(l),  Ex(2),  and  Ex(3).   If  there  are  more  than  three 
resonances,  Ex(3)  through  EX(N)  will  all  be  fitted  with  the 
same  line  width. 

L6=0 
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Ly  Target  geometry: 

=0  Transmission  -  normal  thickness  entered  on  card  #3. 

=1  Reflection  -  normal  thickness  entered  on  card  #3. 

=-1  Effective  thickness  entered  on  card  #3. 

Lg  Output  control: 

=0  Maximum  output . 

=1  Minimal  output  -  no  tables  and  no  intermediate  calculations. 

=-1  Limited  output  -  no  tables. 

Lg  Plotting  control: 

=0  Results  are  not  plotted. 

=2  Results  are  plotted. 

L,„  Repetitions  of  the  NAW  procedure. 

=0  No  repetitions  are  performed. 

>0  In  this  case  it  is  possible  to  evaluate  the  data  with  new 

parameters.   The  new  parameters  will  be  entered  on  additional 
cards  placed  after  the  data. 

L   _L  ,   Inelastic  spectrum  fits  iteration  control.   Only  override 
to  the  default  options  need  be  entered. 

LU=2 

L-^2   Maximum  number  of  iterations;  default  is  L]o=30. 

L,^   Minimum  number  of  iterations;  default  is  L,~=5. 

L14=6   . 

L,,-=0  The  radiation  tail  will  be  calculated  using  Ginsberg's  method. 

<0  The  radiation  tail  will  be  calculated  using  Schiff's  method. 

L,,-L10   Radiation  tail  fit  control.   Only  override  options  need  be 
entered. 

L-j^   Analytic  representation  of  the  theoretically  calculated 

radiation  tail;  default  is  L.  =1. 

lb 

=1  SSW  =  P,  +  P„A  +  P^A2  +  P,A3  +  P.A4 
1     I  J      4      b 

=2  SSW  =  P1  +  P2A  +  P3A2  +  P^E  +  P5E2 
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=3  SSW  =  P2  +  P2A  +  P3E  +  P4E2  +  P5E3 
=4  SSW  =  PjA  +  P2E  +  P3E2  +  P4A2  +  P5E3 
=5  SSW  =  Pjl  +  P2U  +  P3U2  +  P4U3  +  P5U^ 

P-  will  be  determined  internally  by  NAW. 

A  =  l/(Eel  -  E) 

U  =  E  -  E^ 

E,  =  GSSW„  +  (GSSW„  -  GSSW„)/L1Q 

D  U  O  U     lo 

<0   The  radiation  tail  parameters  P.  and  E,  will  be  specified 
on  card  #15. 

L17   The  number  of  parameters  for  the  radiation  tail  fit 
(3<L17<5);  default  is  L]?=4. 

L10    See  Li/-;  default  is  Lio=2. 

L,q>0  Error  ellipsoid  control  values  will  not  be  calculated  near 
the  fit. 

L2Q=0 

L2^>0  The  calculation  of  the  radiation  tail  will  be  modified  by 
the  factors  entered  on  card  #16. 

L  >0  Inelastic  spectrum  results  will  be  punched  out  onto  cards. 

8.  Ex(l),  Ex(2),  ...,  Ex(7),  LE  (7F10.5.I10) 

Resonance  excitation  energies  in  MeV;  enter  those  with 
variable  position  first.   For  more  than  seven  resonances 
set  LEX)  and  then  continue  listing  the  additional  reso- 
nances on  successive  cards. 

8a.   Ex(8)  ...  Ex(14),  LE 

8b.   Ex(15)  ...  Ex(21),  LE 

8c.   E  (22)  ...  E  (25) 

9.  V(l),  V(2),  ...,  V(7),  LBR  (7F10.5.I10) 

Resonance  line  widths.   For  more  than  seven  resonances  set 
LBR>0  and  continue  on  successive  cards  as  on  card  //8.   If 
some  or  all  V(I)  are  zero,  then  the  line  width  used  will  be 
the  dispersion  of  the  spectrometer  counting  system. 

10.  P(l),  P(2),  ...,  P(7),  LP  (7F10.5.I10) 

The  ?£  factors  used  when  L^=2  on  card  #7.   Set  LP>0  for 
more  than  seven  factors  and  continue  as  for  card  #8. 
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11.  Z,  A,  XRAD,  RMS,  R,  G,  C,  T  (8F10.5) 
Z     Atomic  number. 

A     Mass  number. 

XRAD  Radiation  length  in  g/cm^ 

RMS   Root  mean  square  nuclear  radius  in  fm.   Used  only  when  N=2 
on  card  //14,  otherwise  left  blank. 

R,  G  When  N=0  or  N=l  on  card  #14,  R  and  G  are  interpreted  as 

the  parameters  in  the  Helm  expression  for  the  form  factor. 
When  N=3,  RsCRT  and  G=TRT.   When  N=2  or  N=4,  R  and  G  are 
left  blank. 

C,  T  The  fermi  charge  distribution  parameters.   Since  these  val- 
ues are  used  as  the  parameters  for  the  trapezoidal  charge 
distribution  potential  used  in  the  radiation  tail  calcula- 
tions, the  values  must  be  entered  even  when  N/4. 

12.  XMAS,  H  (2F10.5) 
XMAS   Nuclear  mass  equivalent  in  MeV. 

H     Isotope  percentage  (100%  =  1.00). 

13.  A,  B,  C,  D,  E  (5F10.5) 

The  A  . . .  E  are  the  coefficients  for  the  calculation  of  the 
elastic  form  factor  using  a  power  series  expansion  onq  : 

FFel  =  A  +  Bq2  +  Cq4  +  Dq6  +  Eq8 

When  A=0  the  elastic  form  factor  will  be  calculated  using 
the  Born  Approximation. 

.14.   £,  N  (I3.24X.13) 

A     Mass  number. 

N     Selection  of  the  nuclear  model  to  be  used  in  the  Born 
Approximation  to  calculate  the  form  factor: 

=0,1  F  =  3(j1(qR)/qR)exp(-G2q2/2)   Helm. 

=2   F  =  1  -q2  <R2>/6  +  q4<R2>2/120  +  q6<R2>3/5040  Reihe. 
=3   F  =  (1  -  CRT*q2)exp(-  TRT*q2)   Shell  model. 
=4   The  form  factor  will  be  calculated  using  the  Fermi  model. 
15.   PI,  P2,  P3,  P4,  P5,  Eb  (6F10.5) 

This  card  is  read  only  if  L^^<0. 
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P1-P5  Radiation  tail  parameters,  see  L^,-  for  the  expressions. 
Ejj    This  value  is  read  only  when  L^g=-5. 
16.   XSSW,  XM,  XB,  LSSW  (3F10.5.I3) 

This  card  is  read  only  if  Lo^X). 

XSSW  Factor  multiplying  the  area  under  the  elastic  peak. 
XM    Factor  multiplying  the  Moller  term. 
XB    Factor  multiplying  the  Bremsstrahlung  term. 


The  equation  for  the  radiation  tail  is  now  given  by: 
(Max  I) 
F2  a 


SSW  =  (--T-  ■■' '   +  XM(Moller  term)  +  XB(Brehm.  term) ) • (XSSW)AE 


M 
LSSW>0   The  calculated  radiation  tail  will  be  punched  onto  cards. 

The  data  cards  now  follow.   The  elastic  and  inelastic  data  must 
be  separated  by  a  card  having  the  data  energy  step  width  entered 
beginning  in  column  72.   Repeated  evaluations  of  the  same  set  of 
data  (i.e.,  Lj,q>0)  must  have  the  following  set  of  cards  placed 
after  the  data  cards  -  one  set  for  each  repetition  to  be  performed. 

1)  L,   ...  L,q  see  card  #7 

2)  E_.   E  see  card  //5 

o '    u 

3)  E  (1.)  ...  E  (N)       see  card  //8 

4)  V(l)  ...  V(N)         see  card  #9 

5)  P(l)  ...  P(N)         see  card  #10 
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E.   MARPHASE 

MARPHASE  uses  the  Rawitscher-Fischer  phase  shift  code  YALERF  developed 
at  YALE  [Ref.  8J  to  calculate  the  elastic  electron  scattering  cross  sections, 
Mott  cross  sections,  and  form  factors.   A  parameter  set  consists  of  two 
cards. 

1.   Z,  A,  PI,  P2,  P3,  MODEL  (5F10.5,I10) 

Z     Atomic  number. 

A     Mass  number. 

P3=0 

M0DEL=1   Fermi  model  of  the  nucleus  is  assumed  with 


Pir)= 


I 


I  «•  exp(r~c/(t/4.0) 


In  this  case  Pl=c  and  P2=t,  both  in  fermis. 
=2   The  Shell  model  of  the  nucleus  is  assumed  with 


pw=['"(xr]^p[-(i)'j 


In  this  case  Pl=  a    and  P2=A. 

2.  E0,  Tl,  T2,  TH,  R0,  LAG  (5F10.5) 
E0    Electron  beam  incident  energy  in  MeV. 

Tl     Starting  scattering  angle  in  degrees. 

T2    Final  scattering  angle  in  degrees. 

TH    Angular  increment  in  degrees. 

R0    The  fitting  radius  pertaining  to  matching  the  Nuclear 
and  the  Coulomb  wave  functions. 

LAG   The  number  of  phase  shifts. 

3.  Another  card  //2  placed  next  will  cause  the  program  to  use  the 
same  information  for  card  #1  to  compute  the  new  combination  of  incident 
energy  and  scattering  angles.   A  blank  card  is  interpreted  to  mean  that 
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a  computation  is  now  desired  for  a  different  set  of  card  #1  parameters  and 
is  followed  by  a  new  parameter  set,  usw. 
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F .   DUKE 

DUKE  calculates  the  theoretical  cross  sections  for  both  electric  and 
magnetic  giant  resonance  transitions.   Either  the  DWBA  or  the  PWBA  may  be 
selected.   A  parameter  set  consists  of  four  cards.   A  subsequent  parameter 
set  consists  of  only  the  last  three  cards.   A  blank  card  signifies  the  end 
of  the  desired  calculations. 

1.  Z,  AX,  Zl,  El,  E2,  ETRAN,  B0RN,  M0RE  (7F10.5.I10) 
Z     Atomic  number. 

AX  Parameter  c  for  the  Fermi  charge  distribution. 

Zl  Parameter  t  for  the  Fermi  charge  distribution. 

El  Starting  energy  in  MeV. 

E2  Final  energy  in  MeV. 

ETRAN  Excitation  energy  in  MeV.   Set  ETRAN=0  for  a  pure  coulomb 
transition  and  set  ETRAN=E1-E2  if  the  transverse  component 
of  the  cross  section  should  be  considered. 

B0RN=O  Distorted  Wave  Born  Approximation. 

=1  Plane  Wave  Born  Approximation. 
M0RE=O  Electric  transition. 

=1  Magnetic  transition. 

2.  L,  KMAX,  RADVAR,  SKIVAR,  ASG,  ZKMAX,  ETRINT,  LM0D  (2110, 5F10. 5 , 110) 
L     Multipolar ity  of  the  transition. 

KMAX>0  Number  of  partial  waves. 

=0  KMAX  is  calculated  internally:   KMAX  =9+  (El-30)/10. 
RADVAR  Position  of  the  transition  distribution  normalized  to  c. 
SKIVAR  Width  of  the  transition  distribution  normalized  to  t. 
ASG=0   Complete  output  is  printed. 
=1,2,3  Limited  output  is  printed. 
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ZKMAX  This  modifies  KMAX,  as  computed  internally,  to  be 
KMAX  =  KMAX  +  ZKMAX. 

ETRINT  Parameter  used  to  determine  the  integration  steps.  ETRINT  = 
35.0  by  default;  increasing  ETRINT  decreases  the  step  width, 
increases  the  accuracy  and  increases  the  execution  time. 

LM0D=O  Transition  distribution  function  will  be  approximated  by 

the  derivative  of  the  ground  state  fermi  distribution  func- 
tion. 

=1  Transition  distribution  function  will  be  approximated  by  the 
ground  state  fermi  distribution  function. 

3.  HH1,  HH2,  NST1,  NST2 ,  NST3 ,  IST1  (2F10.5 ,4110) 

This  card  is  best  given  as  a  blank  card,  DUKE  then  determines  the 
integration  parameters  using  the  internal  SUBPROGRAM  INDEX;  see 
below  for  assumed  values. 

HH1     Integration  stepwidth  inside  the  potential  in  fermis. 

HH2  Integration  stepwidth  outside  the  potential  in  dimensionless 
units. 

NST1    Length  of  the  internal  region  in  units  of  AX. 

NST2   Length  of  the  external  region  in  dimensionless  units  over 
ETRINT. 

NST3   The  values  of  the  radial  wave  functions  will  be  computed  in 

three  regions  delimited  by:   NST1,  HH1  +  NST2,  HH2  +  i'NST3,  HH2 ; 
where  i  =  1,  2,  . . . ,  IST1. 

When  this  card  is  left  blank  then: 

HH1  =  .25  fm. 

HH2  =  0.3"hc/ y/El-E2' 

NST1  =  3AX/HH1 

NST2  =  (30.  +  (ETRAN/E1)- ETRINT) /HH2 

NST3]  Mean  of  the  integrated  radial  wave  function  within  one  wave- 
IST1J  length  of  the  oscillation  with  the  largest  wavelength.   IST1  =  6, 
and  NST3  =  wavelength/ (HH2- IST1) 

4.  THEMN,  DELTH,  N0,  STANZ ,  NZDG  (2F10. 5 , 1110 ,F10. 5110) 

The  cross  section  for  a  parameter  set  will  be  calculated  for  all 

angles  9  =  THEMN  +  i- DELTH,  i=2,  1,  . . . ,N0.   KMF,  STANZ,  and  KZDG 

are  test  parameters  and  have  no  effect  when  left  blank.  If  this 

card  is  left  blank  then  THEMN  =  30,  DELTH  -  15,  and  N0  =  10. 
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G.   BCALC 

BCALC  is  used  to  calculate  the  transition  strength,  B.   While  written 
specifically  for  use  with  75  MeV  and  90  MeV  incident  energies,  it  is  easily 
changed  to  accommodate  other  values.   BCALC  uses  the  following  data  cards: 
one  card  with  the  number  of  75  MeV  experimental  points;  a  card  for  each  such 

point  containing  the  scattering  angle  in  degrees  and  the  experimental  value 

p 
of  F  ,  FORMAT  (13 ,7X,E10. 4) ;  a  similar  set-up  for  the  90  MeV  experimental 

2 
points;  and  lastly,  the  theoretical  values  for  scattering  angle  and  F  . 

For  each  transition  type  to  be  evaluated  a  set  of  theoretical  values  must 
be  supplied  -  75  MeV  values  followed  by  90  MeV  values.   The  output  from 
BCALC  is  a  listing  of  all  experimental  values  with  the  corresponding  the- 
oretical values  printed  next  to  them,  the  computed  value  for  B,  the  com- 
puted value  of  the  quantity  VAR  (see  Section  IV-A) ,  and  finally  a  listing 

9 

of  the  theoretical  value  of  F  multiplied  by  B. 
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work,  measurements  were  made  to  identify  the  giant  nultipole  resonances  of 

1  Q  "7 

^'Au.   Five  inelastic  scattering  spectra  were  obtained  and  analyzed  t'i 
beam  energies  and  angles  of  75  McV  /  75°,  75  MeV  /  90",  75  MeV  /'  120°, 
90  MeV  /  75°,  and  90  MeV  /  90°.   In  each  spectrum,  giant  multipole  resonances 
were  observed  at  excitation  ei       "f  7.1,  9.1,  10.9,  and  14  MeV.   ] 
strengths  and  multipolar!!  ind  to  1       '  1  '  .'  ""■■•' fin'* . 
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B(E2)=8,000e2fm4 ,  B(E2)=9 ,000e2 fm4 ,  and  B(El)=90e2 fm2 ,  respectively.   The 
full  widths  at  half  maximum  for  these  resonances  were  found  to  be,  respect- 
ively, 3.1,  2.9,  3.1,  and  4.9  MeV. 

The  resonances  found  at  10.9  and  14  MeV  confirm  the  only  other  measure- 
ments made  on  gold  f"4j.   The  resonance  found  at  7.1  MeV  agrees  with  those 
findings  concerning  the  magnetic  character  of  the  transition  but  differs 
in  the  multipolarity  assignment;  however,  our  assignment  may  also  be  con- 
sidered tentative  since  insufficient  data  were  obtained  at  the  backward 
scattering  angles.   The  resonance  found  at  9.1  MeV  has  not  been  previously 
reported. 
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